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A search for an excess of muon-neutrinos from dark matter annihilations in the Sun has been 
performed with the AMANDA-II neutrino telescope using data collected in 812 days of livotime 
between 2001 and 2006 and 149 days of livetime collected with the AMANDA-II and the 40-string 
configuration of IceCube during 2008 and early 2009. No excess over the expected atmospheric 
neutrino background has been observed. We combine these results with the previously published 
IceCube limits obtained with data taken during 2007 to obtain a total livetime of 1065 days. We 
provide an upper limit at 90% confidence level on the annihilation rate of captured neutralinos in the 
Sun, as well as the corresponding muon flux limit at the Earth, both as functions of the neutralino 
mass in the range 50 GeV-5000 GeV. We also derive a limit on the neutralino-proton spin-dependent 
and spin-independent cross section. The limits presented here improve the previous results obtained 
by the collaboration between a factor of two and five, as well as extending the neutralino masses 
probed down to 50 GeV. The spin-dependent cross section limits are the most stringent so far 
for neutralino masses above 200 GeV, and well below direct search results in the mass range from 
50 GeV to 5 TeV. 

PACS numbers: 95.35.+d, 95.30.Cq, 11.30.Pb 



I. INTRODUCTION 

There is an impressive corpus of astrophysical and cos- 
mological observations that indicate that a yet unknown 
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form of non-luminous matter constitutes about 23% of 
the content of the universe. This matter can not be bary- 
onic though, as the abundance of baryons in any form 
is strongly constrained by the inferred primordially syn- 
thesized abundances of deuterium and helium, as well 
as the observed small-scale anisotropy in the cosmic mi- 
crowave background (see, e.g., the review in Rcf. [1]). 
However, the observational limits on non-baryonic dark 
matter in the form of relic stable particles from the Big 
Bang are not so strong and, indeed, a wealth of mod- 
els exist that propose candidates with interaction cross 



sections, masses and relic densities which are compati- 
ble with observations. Generically, any particle physics 
model that provides a stable weakly-interacting massive 
particle (WIMP) is of interest from the point of view of 
the dark matter problem. One class of models that are 
extensively studied for their interest in particle physics 
are supersymmetric extensions of the Standard Model. 
Several flavours of supersymmetry are currently under 
scrutiny at the LHC, and both ATLAS and CMS have 
already probed parts of the parameter space of specific 
benchmark models [2, 3]. In this paper we focus on 
the Minimal Supersymmetric extension to the Standard 
Model (MSSM), which provides a WIMP candidate in 
the lightest neutralino, XiJ ^ linear combination of the su- 
persymmetric partners of the electroweak neutral gauge- 
and Higgs bosons. Assuming _R"parity conservation, the 
neutralino is stable and a good dark matter candidate. 
Accelerator searches and relic density constraints from 
WMAP data allow a lower limit to be set on the mass of 
the MSSM neutralino [4]. Typical lower limits for m^o 
from such studies lie around 20 GeV, depending on the 
values chosen for tan/3, the ratio of the vacuum expecta- 
tion values of the two neutral Higgses. Theoretical argu- 
ments based on the requirement of unitarity set an upper 
limit on m^o of 340 TeV [5] . Within these limits, the al- 
lowed parameter space of minimal supersymmetry can 
be exploited to build realistic models which provide relic 
neutralino densities of cosmological interest to address 
the dark matter problem. 

Relic ncutralinos in the galactic halo can lose energy 
through scatterings while traversing celestial bodies, like 
the Sun, and become gravitationally bound into orbits in- 
side. The build-up of the accreted particles is limited by 
annihilations, which ultimately create high-energy neu- 
trinos [6]. In this paper we present two independent 
searches for a neutrino flux from the annihilations of ncu- 
tralinos captured in the center of the Sun, performed with 
the AMANDA-II and IceCube neutrino telescopes. The 
overwhelming majority of triggers in AMANDA-II and 
IceCube, ©(lO^^/year), are due to atmospheric muons 
reaching the depth of the array. The analyses are there- 
fore based on the search for up-going muon tracks from 
the direction of the Sun when it is below the horizon, and 
have been optimized in order to maximize the sensitivity 
to the predicted signal from MSSM ncutralinos. The re- 
sults can, however, be interpreted generically in terms of 
any other dark matter candidate which would produce a 
similar neutrino spectrum at the detector. The first anal- 
ysis (Analysis A) uses the data taken with AMANDA- 
II during 812 days of livetime in stand-alone operation 
between March 2001 and October 2006. In 2007 the de- 
tector was switched to a new data acquisition system 
(DAQ) which included full waveform recording, and this 
allowed us to integrate the detector as a subsystem of Ice- 
Cube. Analysis B uses 149 days of livetime collected be- 
tween April 2008 and April 2009 (when AMANDA-II was 
decommissioned) using the new AMANDA-II DAQ and 
data from both detectors, AMANDA-II and IceCube. 




FIG. 1. Surface location of the IceCube 40-string configura- 
tion (filled squares) and AMANDA-II (filled circles) 



We have previously published a search for dark matter 
accumulated in the Sun using data taken in 2007 with 
the 22-string configuration of IceCube [7] . In Section VI 
we combine these three independent analyses (Analysis 
A, Analysis B and the results from Ref. [7]) which cover 
the period 2001-2008, presenting competitive limits on 
the muon flux from neutralino annihilations in the Sun 
and limits on the spin-dependent and spin-independent 
neutralino-proton cross section. 



II. THE AMANDA-II AND ICECUBE 
DETECTORS 

The AMANDA-II detector consisted of an array of 677 
optical modules (OM) deployed on 19 vertical strings at 
depths between 1500 m and 2000 m in the South Pole 
ice cap. The optical modules consisted of a 20 cm diam- 
eter photomultiplier tube housed in a glass pressure ves- 
sel connected through a cable to the surface electronics, 
where the photomultiplier pulses were amplified, time- 
stamped and fed into the trigger logic. The inner ten 
strings used electrical analog signal transmission, while 
the outer nine strings used optical fiber transmission. 
The strings were arranged in three approximately con- 
centric circles of 60 m, 120 m and 200 m diameter respec- 
tively. The vertical separation of the optical modules in 
strings one to four was 20 m, and 10 m in strings five 
to nineteen. Several calibration devices (a 337 nm N2 
laser and three DC halogen lamps, one broadband and 
two with filters for 350 and 380 nm) were deployed at 
several locations in the array. Additionally, a YAG laser 
calibration system was set up on the surface, able to send 
light pulses through optical fibers to a diffuser located at 
each optical module. 

The original AMANDA-II data acquisition system 
could only measure the maximum amplitude of the pho- 
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tomultiplier pulses and only eight hits could be buffered 
at a time. From 2003, each channel was connected to 
a Transient Waveform Recorder (TWR), a flash ADC 
that samples at 100 MHz with 12 bit resolution, being 
able to collect the full photomultiplier waveform in each 
channel [8]. The TWR system was run in parallel with 
the original DAQ until 2007, when the original DAQ was 
switched off. From this point, until its decommissioning 
in May 2009, AMANDA-II was run only with the TWR 
DAQ system. 

The construction of the IceCube detector began in 
2005 and 40 strings surrounding the AMANDA II array 
had been deployed by 2008. With its denser string con- 
figuration, AMANDA-II played the role of a low-energy 
subdetector of IceCube, and both detectors where run 
in an 'OR' configuration until AMANDA-II was decom- 
missioned. Figure 1 shows the siuface location of the 
AMANDA-II and IceCube strings during the 2008 data 
taking period. The IceCube strings consist of 60 Digital 
Optical Modules (DOMs) separated vertically by 17 m. 
The inter-string separation is 125 m. Contrary to the 
AMANDA-II optical modules, the IceCube DOMs dig- 
itize the photomultiplier signals in situ. They are also 
self-calibrating units, frequently exchanging timing cal- 
ibration signals with the surface electronics. A descrip- 
tion of the IceCube DOM and the DAQ can be found in 
Refs. [9, 10]. 

Both IceCube and AMANDA-II used triggers that se- 
lected events based on the number of hit optical modules 
in a certain time window (at least 24 modules hit within 
2.5 fis in AMANDA-II, and at least 8 modules hit within 
5 /is in IceCube). Additionally, a trigger tailored for low 
energy events required N hit modules in a given time 
window out of M consecutive modules in the same string 
(N/M condition). This trigger was set to 6/9 in strings 
one to four and 7/11 in strings five to nineteen for most 
of the AMANDA II livctime, with a time window of 2.5 
/iS. The condition was 5/7 in all IceCube strings, with a 
time window of 1.5 jis. 

Muons from charged-current neutrino interactions 
near the array arc detected by the Chercnkov light they 
produce when traversing the ice. The hit times along 
with the known detector geometry and the optical prop- 
erties of the ice allow the reconstruction of the tracks 
passing through the detector. A detailed description of 
the reconstruction techniques used in AMANDA and Ice- 
Cube is given in Ref. [11] 



III. SIGNAL AND BACKGROUND 
SIMULATIONS 

The simulation of the neutralino-induced neutrino sig- 
nal was performed using the WimpSim program [12] for a 
sample of neutralino masses (50, 100, 250, 500, 1000, 
3000 and 5000 GeV). Two extreme annihilation chan- 
nels were considered in each case, a soft neutrino chan- 
nel, XiXi ~^ bb and a hard neutrino channel, XiXi ~^ 



W'^W~ , (xiXi T~^T~ was chosen for 50 GcV ncu- 
tralinos, since they are assumed to annihilate at rest and 
this mass is therefore below the W production thresh- 
old). This choice covers the range of neutrino energies 
that would be detectable with AMANDA-II/IceCube for 
typical MSSM models. Note, however, that neutrinos 
with energies above about a few hundred GeV will inter- 
act in the Sun, and do not escape; only the lower energy 
neutrinos from the decays of the products of such interac- 
tions will get out of the dense solar interior. This cutoff, 
rather than the WIMP mass, sets the neutrino spectrum 
for WIMP masses above 1 TeV. The simulated angular 
range was restricted to zenith angles between 90° (hori- 
zontal) and 113°, with the generated number of events as 
a function of angle weighted by the time the Sun spends 
at each declination. WimpSim propagates the neutrinos 
taking into accoimt energy losses and oscillations in their 
way out of the Sun as well as vacuum oscillations to the 
Earth, giving the expected neutrino flux at the location 
of the detector for each neutralino mass simulated. The 
neutrino-nucleon interactions in the ice around the de- 
tector producing the detectable muon flux were simulated 
with nusigma [12] using the CTEQ6 [13] parametrization 
of the nuclcon structure functions. 

The background for the dark matter searches 
arises from up-going atmospheric neutrinos and mis- 
reconstructed downward going atmospheric muons. Wc 
have simulated the atmospheric neutrino flux according 
to Ref. [14], using the ANIS program [15], with energies 
between 10 GeV and 325 TeV and zenith angles between 
80° and 180° (vertically up-going). The simulation in- 
cludes neutrino propagation through the Earth, taking 
into account the Earth density profile [16], neutrino ab- 
sorption and neutral current scattering. The simulation 
of atmospheric muons was based on the CORSIKA air 
shower generator [17] using the South Pole atmosphere 
parameters and the Horandel parametrization of the cos- 
mic ray composition [18]. We have simulated 10^"'^ in- 
teractions, distributed uniformly with zenith angles be- 
tween 0° and 90°, and with primary energies, Ep, be- 
tween 600 GcV and 10" GeV. We note that the back- 
ground simulations were not used in the evaluation of 
the actual background remaining at final cut level in the 
analyses presented below, but off-source data was used 
to that end. Background simulations were used as a con- 
sistency check at the different steps of the analyses. 

Muons were propagated from the production point 
to the detector taking into account energy losses by 
bremsstrahlung, pair production, photo-nuclear interac- 
tions and (5-ray production as implemented in the code 
MMC [19]. The Cherenkov light produced by the muon 
tracks and secondaries was propagated to the optical 
modules taking into account photon scattering and ab- 
sorption according to the measured optical properties of 
the ice at the detector [20], as well as the measured dust 
layer structure of the ice at the South Pole [21]. The 
Photonics [22] program was used to this end. 
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TABLE I. ANALYSIS A: For each ncutrahno mass and annihilation channel the table shows: The median of the space angle 
distribution, 'I', the effective volume, Voft, the 90% CL upper limit on the expected signal, /if", and the corresponding 90% CL 
limits on the annihilation rate at the center of the Sun, Fa, on the rrmon flux at the Earth, $p, and on the spin-dependent 
and spin-independent neutralino-proton cross sections, rf^ and a^p. The limits include systematic uncertainties. 



IV. DATA ANALYSIS 

Below we describe the two analyses performed with 
AMANDA-II data in stand-alone mode and with 
AMANDA-II run in coincidence with the 40-string con- 
figuration of IceCube. In both cases the cuts were opti- 
mized on data when the Sun was above the horizon, so 
the data analysis was kept blind to the actual direction 
of the Sun. 



A. ANALYSIS A 
(AMANDA-only analysis) 

The data set used in this analysis corresponds to a to- 
tal of 812 days of effective livetime, and comprises a total 
of 7.25x10^ events collected when the Sun is below the 
horizon between the beginning of March 2001 and the end 
of October 2006. Data from periods where the detector 
showed unstable behaviour or periods where test or cali- 
bration runs were performed, were removed from the final 
data sample. The remaining events were cleaned of hits 
induced by electronic cross-talk, dark noise or unstable 
modules, and the data set was retriggered to make sure 
that the required amount of physical hits participated in 
building the trigger. The events were first reconstructed 
with two fast first-guess track finding algorithms. Di- 
rect Walk [11] and JAMS [23]. These reconstructions arc 
aimed at identifying muon tracks and give a rough first 
estimate of their direction, using it as a first angular cut 
in order to reduce the data sample, still dominated by 
down-going atmospheric muons at this level. Events with 
zenith angles smaller than 70° as reconstructed by the 
DirectWalk or the JAMS reconstructions were rejected 
at this stage. Two maximum likelihood reconstructions 



were applied to the remaining events. One uses Direct- 
Walk and JAMS as seeds, performing an iterative maxi- 
mization of the probability of observing the actual event 
geometry (hit times and positions) with respect to a given 
track direction. In reality, the negative of the logarithm 
of the likelihood is minimized in order to find the best-fit 
zenith and azimuth. In order to exploit the fact that most 
of the events that trigger the detector are down going 
atmospheric muons, an iterative Bayesian reconstruction 
incorporating the known atmospheric muon zenith angle 
distribution as prior, was also performed. A comparison 
between the likelihoods of the standard and the Bayesian 
reconstructions can then be used to evaluate the likeli- 
hood that an event is down-going or up-going. 

The reconstructed events were then processed through 
a series of more stringent angular cuts on the direction of 
the tracks as obtained with the likelihood fits (zenith an- 
gle > 80°), to reduce the atmospheric muon background 
and retain as much of a potential signal as possible. At 
this level, only a fraction of 3x 10~^ of the data and the 
simulated atmospheric muon background survive, while 
69% of the simulated atmospheric neutrino background 
was kept. Between 36% and 78% of the ncutralino signal 
survived these cuts, depending on the neutralino mass 
and annihilation channel. 

The final event selection step was performed using a 
boosted decision tree (BDT) classifier [24]. We have se- 
lected 21 variables that showed good separation power 
and a correlation below 65% between any pair. The vari- 
ables used for tlic^ classification scheme can be grouped in 
two major classes: variables related to the hit topology of 
the event and variables related to the quality of the track 
reconstruction. Among the first class there arc variables 
like the number of hit optical modules, the number of 
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non-scattcrcd, or 'direct', hits"'^, the number of strings 
with hits within a 50 m radius cyhnder around the track, 
and the number of strings with direct hits in the same 
cyhnder, the center of gravity (c.o.g.) of the spatial po- 
sition hits, the distance of the c.o.g. to the geometrical 
center of the detector, the length of the direct hits pro- 
jected onto the direction of the track, the smoothness of 
the distribution of direct hits along the track and vari- 
ables related to the probability of detecting a photon in 
the hit modules given the reconstructed track hypoth- 
esis. Among the variables related to the reconstruction 
quality, we have used the angle of the standard maximum 
likelihood reconstruction, the difference of log-likelihoods 
between the standard reconstruction and the Bayesian re- 
construction and a measure of the angular resolution of 
the first-guess reconstructions. 

In order to exploit the differences in the final muon en- 
ergy spectra at the detector produced by the annihilation 
of neutralinos of different mass, we trained BDTs sepa- 
rately for each ncutralino mass and annihilation channel, 
using the same 21 variables in each case. The signal 
training samples consisted of 50% of each of the simu- 
lated signal samples. Data when the Sun is above the 
horizon were used as the background sample. The BDT 
classifies the given data as background or signal accord- 
ing to a continuous parameter that takes values between 
1 (signal-like) and -1 (background-like). A cut on the 
BDT output was chosen for each model as to maximize 
the discovery potential. Details of the data analysis and 
a complete list of the variables used can be found in 
Ref. [25]. 

After the BDT cut the data sample is reduced by a fac- 
tor between lxlO~^ and 3xl0~^ with respect to trigger 
level, depending on the signal model used for the opti- 
mization, and its purity is closer to the irreducible atmo- 
spheric neutrino background. Wc note however that the 
approach taken in this and the analysis described in Sec- 
tion IV B does not require to reach a pure atmospheric 
neutrino sample, since we use the shape of the normal- 
ized space angular distributions with respect to the Sun 
of both signal and background to build our hypothesis 
testing. Details are given in Section V below. 

1. Systematic uncertainties 

The eff'ect of different systematic uncertainties on the 

signal expectation was evaluated by varying the relevant 
parameters in the Monte Carlo and processing this new 
sample through the same analysis chain as the nominal 
sample. Systematic uncertainties affect the sensitivity of 
the detector to a given signal, a quantity that is char- 
acterized by the efii'ective volume. The effective volume 



^ Hits with a time stamp that corresponds to the time that light 
taJces to travel directly from the track hypothesis to the optical 
module. 



is defined as Vcs = ('T^finai/%on) x Vgcn and calculated 
by generating a given number of signal events rigen in 
a geometrical volume Vgcn around the detector, rifinai 
is the remaining number of events after the analysis cuts 
have been applied to the generated event sample. The ef- 
fective volume depends on neutrino energy, but in what 
follows wc will show the integrated volume over the neu- 
trino spectrum produced by a given neutralino mass. The 
relative uncertainty in the effective volume is then given 

by, 

AKff ^ Vsys - VcS . s 

where Vcs is the effective volume of the baseline analy- 
sis and Vsys tfie effective volume calculated with a given 
assumption for systematic effects. Since we are not rely- 
ing on atmospheric muon and neutrino Monte Carlo to 
estimate the background, we evaluate the effect of sys- 
tcmatics only on the signal expectation. 

Systematic uncertainties can be classified in different 
categories. There arc systcmatics induced from the un- 
certainties in quantities used in the signal Monte Carlo; 
neutrino cross sections, oscillation parameters or muon 
energy losses. The uncertainties in the oscillation pa- 
rameters used to calculate the expected neutrino fiux 
at the detector were taken from Ref. [26] and lead to 
an uncertainty of less than ±3%. Further, uncertain- 
ties from the neutrino-nucleon cross section calculation 
within WimpSim and the simulation of muon energy loss in 
the ice have been estimated to be ±7% [12] and ±1% [19] 
respectively. 

An additional source of systematics is the implemen- 
tation of the optical properties of the ice in the detector 
response simulation, as well as the in-situ optical module 
sensitivity. The AMANDA-II calibration light sources 
were used to measure photon arrival time distributions as 
a function of the relative emitter-receiver distance. Such 
measurements allow us, in principle, to extract an effec- 
tive scattering length and absorption length that charac- 
terize the deep ice where the detector is located. How- 
ever, the ice at the South Pole presents a layered struc- 
ture, with slightly different optical properties due to the 
presence of different concentrations of dust at different 
depths [27]. Moreover the process of melting the ice to 
deploy the strings and the subsequent refreeze of the wa- 
ter column changes the local optical properties of the 
ice in the drill hole. The inverse problem of extracting 
the ice properties from the measured photon arrival time 
information is then a quite difficult one, and different im- 
proved implementations have been applied in AMANDA 
with time. We have used the two most recent ice models 
as an estimation of the uncertainties introduced by the 
different implementations of the ice properties as a func- 
tion of depth. The relative uncertainty induced in the 
detector effective volume by this effect ranges between 
3% and 30%, depending on neutrino energy. 

The uncertainty on the total sensitivity of the deployed 
optical modules (glass plus PMT) also contributes to 
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FIG. 2. Monte Carlo neutrino energy distribution at final 
level of Analysis B, of the events triggering only AMANDA- 
II (full line), only IceCube (dashed line) or both detectors 
(dotted line). 



the uncertainty in the effective volume. Two additional 
Monte Carlo samples with the light collection efhciency 
of each module globally shifted by ±10% with respect 
to the baseline simulation were produced and the events 
passed through the complete the analysis chain. The ef- 
fect on the effective volume lies between 20% and 40%. 
The overall systematic uncertainties in the detector ef- 
fective volume lie between 20% and about 50% depend- 
ing on the neutralino mass model being tested. The 
total uncertainty in VeS has been calculated under the 
assumption that all uncertainties are uncorrelated, i.e., 

Av/v = ^TJ^vmi 



B. ANALYSIS B 
(AMANDA-IceCube combined analysis) 

The data set used in this analysis comprises a total 
of 1.7x10^° events collected between the middle of April 
and end of September 2008, plus the period March 20 
- April 1 2009. Only runs where both AMANDA-II 
and IceCube were active were used, which corresponds 
to a total livetime of 149 days when the Sun was be- 
low the horizon. The surface geometry of AMANDA- 
II/IceCube-40 used in the analysis is shown in Fig- 
ure 1. The denser AMANDA-II array plays the role of 
a low energy detector, as can be seen in Figure 2 which 
shows the energy distribution of neutrinos triggering the 
AMANDA-II array and the IceCube array, at final cut 
level. We have therefore simplified the analysis with re- 
spect to the approach taken in analysis A, and optimized 
the cuts for two energy regions, low neutrino energies and 
high neutrino energies. We used 100 GeV neutralinos an- 
nihilating into bb as a benchmark for the low-energy op- 
timization, and 1000 GeV neutralinos annihilating into 



W^W~ for the high-energy optimization. The opti- 
mization that was finally used on each WIMP model was 
decided at the end of the analysis, based on which one 
achieved the best sensitivity for the given model. 



1. Data analysis 

Data analysis proceeded on slightly different lines for 
the low-energy and high-energy streams. Events that 
had hits in the AMANDA-II detector were reconstructed 
with the JAMS first-guess reconstruction and then with 
an iterative likelihood fit. For events with hits only in the 
IceCube strings, a simple line-fit [11] proved to work well 
and it was used as the first-guess. A series of straight 
cuts on event quality and track direction were performed 
in order to reduce the amount of data while keeping as 
much of the signal as possible, before using a multivari- 
ate classifier for the final separation. The cut variables 
used at this stage were slightly different between the low- 
energy and the high-energy streams, reflecting the differ- 
ent geometry and size of the AMANDA-II and IceCube 
detectors, as well as the fact that we want to use Ice- 
Cube as a veto for the low-energy stream. Events classi- 
fied as belonging to the low-energy stream were required 
to have an AMANDA-II trigger, a reconstructed zenith 
angle larger than 90°, at least 25% of direct hits (with 
a minimum of four direct hits), a distance between the 
first and last direct hit of at least 25 m, and less than five 
hits in any of the IceCube strings. Events in the high- 
energy stream were also required to have a reconstructed 
zenith angle larger than 90° and smaller than 120°, a 
reduced log-likelihood smaller than 6.5, a difference be- 
tween the angles of the likelihood reconstruction and the 
first-guess fit smaller than 40°, at least 3 direct hits, at 
least 2 strings with direct hits and a distance between 
the first and last direct hit of at least 141 m. These cuts 
reduced the data and atmospheric muon simulation by a 
factor of about 2000 in each case, while just reducing the 
low-energy and high-energy signal streams by a factor 
of 9 and 2.4, respectively. 

After these straight cuts, a support vector machine 
(SVM) [28] was used for the final classification of the 
events. We trained two SVMs independently, one for the 
low-energy sample and another one for the high-energy 
sample. The signal training samples consisted of 50% of 
the simulated signal samples considered in each stream. 
Data when the Sun is above the horizon were used as the 
background sample. The choice of variables to use in the 
SVM was done iteratively: from an original group of 35 
variables, the SVM was trained on all variables but one 
at a time. A variable for which the SVM performance 
did not worsen when removed, was discarded. A set of 
12 variables were identified by this method as useful for 
signal and background discrimination in the low-energy 
stream, and ten for the high-energy stream. As in Anal- 
ysis A, the variables are related to the event topology 
and the quality of the track reconstruction. Three of the 
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FIG. 3. Figure of merit of the model discovery potential (upper panels) and sensitivity (lower panels) as a function of SVM 
output of Analysis B. The left plots correspond to the low-energy analysis optimization and the right plots to the high-energy 
optimization. 



variablcis wciic common to both streams: the zenith differ- 
ence between the reconstructed track and the position of 
the Sun, the angular difference between the reconstructed 
track and a line joining the centers of gravity of the first 
and the last (in time) quartiles of the hits, and the hor- 
izontal distance from the detector center to the center 
of gravity of the hits. Among the other variables used 
in a particular stream are the number of hit strings, the 
angular resolution of the track fit and the smoothnciss of 
the distribution of direct hits projected along the track 
direction. Details of the data analysis and the complete 
list of the variables used can be found in Rcf . [29] . The 
SVM classifies the given data as background or signal ac- 
cording to a continuous parameter, Q, that takes values 
between (signal-like) and 1 (background like). A cut 
value of Q = 0.1 was used for both data sets, and it is jus- 
tified in Section V below. A fraction of 2.2x10"^ of the 
data in the low-energy analysis and 4.1 x 10~^ of the data 
in the high-energy analysis remain after this cut, while 
24% and 39% of the signal are retained, respectively. No 
further cuts were applied after the SVM classification. 



2. Systematic uncertainties 

The uncertainties in analysis B were evaluated as in 
the previous analysis, by varying the optical module sen- 
sitivity and the ice model in the Monte Carlo and evalu- 
ating the effect on the effective volume after all analysis 
steps. The values for the uncertainties pertaining to theo- 
retical inputs (neutrino-nucleon cross section, oscillation 
parameters and muon propagation in ice) were taken as 
for analysis A. 

With the deployment of IceCube, the ice properties at 
deeper depths than the AMANDA location needed to be 
evaluated, and a new ice model was developed that cov- 
ered the whole IceCube volume and a better description 



of the layered structure of the ice. Two approaches were 
taken in updating the AMANDA ice properties model- 
ing: extrapolating the AMANDA model based on new ice 
core data and using IceCube flasher data. The difference 
in the effective volume induced by these two approaches 
was taken as an estimation of the uncertainty in the anal- 
ysis due to the ice optical properties. This tmcertainty 
ranges from 2% to 4%, depending on the analysis stream, 
while the DOM sensitivity uncertainty lies between 12% 
and 24%. The uncertainties were added in quadrature to 
obtain the total uncertainty of each of the two streams 
of the analysis. 

V. RESULTS 

The BDT classifiers of Analysis A were trained to 

maximize the model discovery potential (MDP) as de- 
fined in Ref. [30], which measures the possible contri- 
bution of a signal when a given number of background 
events, n^, are expected. Reference [30] gives an analyt- 
ical parametrization for the sensitivity of the experiment 
for a chosen discovery level at a given confidence level. 



f + ^+a^+|V^2 + 4a^ + 4nB 

(2) 

where represents the signal efficiency. This figure of 
merit was chosen to represent the strength of the signal 
flux needed for a 5a signiflcance discovery (a=5) at 90% 
confldence level (6=1.28). 

The SVMs in Analysis B were trained to optimize the 
sensitivity (the median upper limit of the number of sig- 
nal events at 90% confidence level, /Ig" divided by the de- 
tector effective volume, Ves) to the signal, although such 
optimization was found to be equivalent to maximizing 
the discovery potential. This is illustrated in Figure 3, 
where lower panels show the sensitivity, while the upper 
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TABLE II. ANALYSIS B: For each neutralino mass and annihilation channel the table shows: The median of the space— angle 

distribution, ^, the effective volume, Voh, the 90% CL upper limit on the expected signal, /if'', and the corresponding 90% CL 
limits on the annihilation rate at the center of the Sun, Fa, on the muon flux at the Earth, <1>m, and on the spin-dependent 
and spin-independent neutralino-proton cross sections, cr^ and cr|p. The limits include systematic uncertainties. 



panels show the figure of merit of the discovery poten- 
tial, both as a function of the cut on the SVM output Q. 
The left plot corresponds to the low-energy optimization 
and the right plot to the high-energy optimization. The 
figure shows that the chosen cut on the SVM {Q = 0.1), 
based on the sensitivity plots, lies very close to the max- 
imum of the discovery potential curves. 

We use the shape of the space angle distribution of 
the final data samples with respect to the Sun to build 
a confidence interval for n^, the number of signal events 
compatible with the observed data distribution at a given 
confidence level, that we will choose to be 90%. The plots 
in Figure 4 show an example of the angular distribution 
of the remaining data after all cuts, compared with the 
angular distribution of the expected background obtained 
from time-scrambled data from the same declination of 
the Sun. Angles shown are space angles with respect to 
the Sun position, and the different plots correspond to 
Analysis A and Analysis B as indicated in the caption. 

A probability density function for the space angle 
can be constructed as follows, 

/(^IMs) = + f 1 - — ) ./b(V) (3) 

^obs \ ^obs / 

representing the probability of observing a given angle 4' 
when /ig signal c!vc;nts are present among the total num- 
ber of observed events riobs- The functions fsitp) and 
fBi'fp) are the angular probability density functions of 
signal and background respectively, obtained by fitting 
and normalizing the corresponding angular distributions. 
The likelihood of the presence of /ig signal events in an 
experiment that observed exactly riobs events, can then 
be expressed as follows 

"obs 

= n /(V'ilMs), (4) 

i=l 



To define confidence intervals based on Eq. (3) we use 
the likelihood-ratio test statistic 7^(/Xs) = £(/Us)/£(/is) 
as proposed by Feldman and Cousins [31], where /Is is 
the result of best fit to the observed ensemble of space 
angles. For each ^s, ^{f^s) < ^(J^s) and TZ{fis) < 1- 
Through a series of pscudo experiments with varying val- 
ues of a critical value TZ^°{iJ.s) can be found such as 
lnR.{fis) >ln7?.^''(/Xs) for 90% of the cases. A confidence 
interval [/xiow, Mup] at 90% confidence level can then be 
calculated as [/Uiow; Atup]={A* I 1ii'7^(Ms) >1ii7?.^°(/Us)}- 

In the absence of a signal, /lup is the 90% confidence 
level limit on the number of signal events, /i^'^. This 
limit can be directly transformed into a limit on the vol- 
umetric rate of neutrino interactions in the detector due 
to a signal flux, r,^^^, since we would expect a num- 
ber Us = ^ly^fiVes ■ tiive of neutrinos in the livetime tuve- 
Therefore 

90 

VeS ■ Hive 

This volumetric interaction rate is directly propor- 
tional to the neutralino annihilation rate in the Sun, Fa, 
through 

^"^1^ = ^^^2 ^^'^ i^U- ^ ^^^^ I ^'') 

^Posc(m,^)E^k(^) (6) 

where Dq is the distance to the Sun, a^M the neutrino- 
nucleon cross section (above a given muon energy thresh- 
old, E^'^, taken as 10 GeV), pN the nucleon number den- 
sity of the detector medium, Posc{^^,i) the probability 
that a produced neutrino of flavour i oscillates to flavour 
H before reaching the detector (including the probability 
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FIG. 4. The space angular distribution with respect of the position of the Sun (first 15 degrees) of the remaining data events 

in two analyses, after all cuts (dots). The histogram in all plots represents the expected background distribution obtained from 
time-scrambled data in the same declination of the Sun, with the la Poisson uncertainty shown as the shaded area. Upper 
row: Analysis A: The left plot shows the case of optimizing the analysis for 100 GcV neutralinos and bb annihilation channel 
as an example of a low-energy optimization. The right plot shows the case of 1000 GeV neutralinos and W^W' annihilation 
channel as an example of a high-energy optimization. Lower row: Analysis B: The left plot shows the case of optimizing the 
analysis for 100 GeV neutralinos and bb annihilation channel (the low energy optimization). The right plot shows the case of 
1000 GeV neutralinos and W^W' annihilation channel (the high energy optimization) 



that the neutrino escapes the dense solar interior), Bk is 
the branching ratio for annihilation into channel K, and 
dN^ I dE^ the number of neutrinos of flavour i produced 
per annihilation and unit of energy from channel K. The 
branching ratios to a given channel are the only unknown 
quantities in the above equation. They depend on several 
unknown SUSY parameters, i.e. the composition and the 
mass of the neutralino. To be able to make concrete pre- 
dictions and simplify the way the results are presented, 
we have performed the analysis on two annihilation chan- 
nels per simulated mass, W^W~ (or t+t~ if to^ < mw) 
and bb, assuming 100% branching ratio to each channel 
in turn. For a given mass, these channels produce the 
hardest and softest neutrino energy spectra respectively 
and are taken as representative of the range of possible 
outcomes if nature chose the MSSM neutralino as dark 
matter. 



Equivalently to Eq. (6), we can relate the annihilation 

rate and the muon flux at the detector, above the 
energy threshold E^^'^ as 

daUE,,El) /diVKX 

e^osc(m,oi:^k(^^j (7) 

where da^{Ei,, E^)/dE^ is the differential neutrino cross 
section for producing a muon with energy E^^ from a neu- 
trino of energy E,^. The term X{E^) denotes the range 

of the muon produced with energy E^, and takes into 

account energy losses between the production point and 
the detector location. For each mass and annihilation 
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TABLE III. Results from the combination of Analysis A, Analysis B and the results of Rcf. [7], covering 1065 days of livetirnc 
between 2001 and early 2009. For each neutralino mass and anuiliilation channel the table shows: The 90% CL upper limit on 
the expected signal, /ig", and the corresponding 90% CL limits on the armihilation rate at the center of the Sun, Fa, on the 
muon flux at the Earth, <I>^, and on the spin-dependent and spin-independent neutralino-proton cross sections, a^p and cr^p. 
The limits include systematic uncertainties. 



channel considered, the integrals have been performed 
with WimpSim, providing a relationship between the ex- 
perimentally obtained limit on F^^^ and derived limits 
on Fa and 

If we assume equilibrium between the capture and 
annihilation rates, the annihilation rate is proportional 
to the neutralino-proton scattering cross section, which 
drives the capture. Under the further assumption that 
the capture rate is fully dominated either by the spin- 
dependent (SD) or spin-independent (SI) scattering, we 
can extract conservative Umits on either the SD or the 
SI neutralino-proton cross section from the limit on Fa- 
We have followed the method described in Ref. [32] in 
order to extract limits on the SD and SI cross sections. 
This conversion introduces an additional systematic un- 
certainty in the calculation of the cross sections, due to 
the element composition of the Sun, effect of planets on 
the capture of halo WIMPs and nuclear form factors used 
in the capture calculations. These effects influence the 
SI and SD calculations differently, and are discussed in 
Ref. [32] . The uncertainties introduced by the conversion 
are small for the spin-dependent cross section, 2%, and 
larger for the spin-independent cross section, 25%. We 
have added these contributions in quadrature in order to 
obtain the total systematic uncertainty on these quanti- 
ties that we use in this work. We just note that recent 
studies favour a somewhat higher value of the dark mat- 
ter density, po, closer to 0.4 GeV/cm^ [33], than the stan- 
dard value of 0.3 GeV/cm^ assumed in our conversion. 
Since variations in the dark matter density translate in- 
versely into the calculated WIMP-proton cross section, 
the conservative assumption on po makes our limits con- 
servative as well. A recent discussion on the effects of 
uncertainties on the structure of the dark matter halo. 



the dark matter velocity dispersion and the effect of the 
gravitational influence of planets on the capture rate of 
local dark matter has been presented in detail in Ref. [34] . 

The results on the quantities described above for each 
neutralino mass and annihilation channel considered are 
presented in Table I for Analysis A and in Table II for 
Analysis B. The tables show the median of the space 
angle distribution with respect to the Sun, the effective 
volume for each analysis optimization, the 90% CL limits 
on the expected signal, on the annihilation rate at the 
center of the Sun, on the muon flux at the Earth, and 
on the spin-dependent and spin-independent neutralino- 
proton cross sections. 



VI. COMBINATION OF RESULTS 

Given that the data samples used in the analyses pre- 
sented in this paper are independent, we can combine the 
results in a statistically sound way. Since IceCube has al- 
ready published an analysis using the same method on 
another independent data set, the data collected with the 
22-string detector in 2007, we can combine these three 
analyses and cover a total livetime of 1065 days, stretch- 
ing from March 2001 to April 2009. We use the combined 
likelihood constructed from the likelihoods of the three 
analyses, 

£(Ms) = 'Ci(/XsWi)£2(/^sW2)'C3(/XsW3) = 
nf"' n°^' nf" 

Yl /l(V'i,l|/"sWl)x]J /2(V'i,2|MsW2)x JJ /aC'^i'i.slMsWs) 
i=l i=l i=l 

(8) 
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where fj,s is now weighted by the Hvetime iuvci and ef- 
fective volume VcS^i of each analysis through the weights 
uji, defined as 



0.05 <n/i' <0.20 



MSSM muon flux limits 



,z ^livc,z 



(9) 



We have chosen to show conservative limits and there- 
fore we have used in Eq. 9 the effective volumes of each 
analysis reduced by its la systematic uncertainty. A 90% 
confidence limit for fig is then obtained using the same 
procedure as for the single-analysis case explained in Sec- 
tion V. The combined limit on F^-i.^ is now given by 



,90 



1 ^cff j' ^livc,j 



(10) 



The calculation of the combined limits on the anni- 
hilation rate and muon flux follows from T^^^ as in 
Eqs. (6) and (7). The results of this procedure are 
shown in Table III and in Figures 5 and 6. The fig- 
ures show the 90% CL limits on the muon flux, on the 
spin-dependent and on the spin-independent neutralino- 
proton cross sections, compared with current limits from 
other experiments. The shaded area shows the allowed 
parameter space of the 7-parameter MSSM, obtained by 
a grid scan using DarkSusy [35]. In order to choose 
only allowed models we have taken into account cur- 
rent experimental limits on the neutralino mass from 
LEP [1] and limits on the WIMP cross section from 
the CDMS [36] and XENON [37] direct detection exper- 
iments. We have allowed for a generous range of values 
of the dark matter relic density D.-^h? around the favored 
value of WMAP [38] , accepting models in the scan which 
predicted values of Q,^h? between 0.05 and 0.2. 

An independent analysis using the point -source search 
techniques described in Ref. [39] has been performed, us- 
ing the Sun as another point source, and has been pre- 
sented in Ref. [40] . The analysis used the cuts developed 
for the point-source search without any further optimiza- 
tion. The only difference being that the estimated energy 
of the events was not included in the likelihood test in 
order to avoid the bias to the E~'^ spectrum used in the 
optimization of the point-source analysis. Since this data 
set (2000-2006) practically overlaps with that of Analysis 
A, we have not included the results of such an approach 
in the combination presented above, but the results of 
the point-source analysis provide a useful confirmation 
of the robustness of the limits presented in this paper. 



VII. DISCUSSION 

We have presented two independent analyses search- 
ing for neutralino dark matter accumulated in the cen- 
ter of the Sun using the AMANDA-II and the 40-string 
IceCube detectors, covering different data taking peri- 
ods. We have combined the obtained 90% confidence 



IC/AMANDA 2001-2008, 44 
IC/AMANDA 2001-2008, ir+ W 



r7s,<4" CDMS(2010)-I-XENON100(2011) 
IC86 180 days sensitivity, 
Super-K2011, 44 
Super-K 2011, B'+ ir 




10' 10' 
WIMP mass (GeV) 



FIG. 5. 90% confidence level upper limits on the muon flux 
from the Sun, from neutralino annihilations as a function 
of neutralino mass. The results from the analyses presented 
in this paper are shown as the black dots joined with lines to 
guide the eye (solid and dashed for the W~^W~ and the bb an- 
nihilation channels respectively). The shaded area represents 
the allowed MSSM parameter space taking into account cur- 
rent accelerator, cosmological and direct dark matter search 
constrains. The red curve shows the expected sensitivity of 
the completed IceCube detector. Super-K results [42] are also 
shown for comparison. 



level limits on the muon flux and neutralino-proton spin- 
dependent and spin-indepenent cross section with the 
previously published limits from the 22-string IceCube 
configuration, to obtain limits corresponding to a total 
livetime of 1065 days. We compare these results to the 
currently allowed MSSM~7 parameter space, obtained by 
a grid scan using DarkSusy. The limits on the muon flux 
are the most constraining so far for neutralino masses 
"^x" ^ GeV. Assuming that the neutralinos consti- 
tute the bulk of dark matter in the Galaxy and that they 
annihilate producing a hard neutrino spectrum at the 
detector, the combined results of AMANDA-II and Ice- 
Cube presented in this paper start to probe the allowed 
7-parameter MSSM space. This comparison with the 
MSSM is, by construction of these analyses, limited to 
a model-by-model rejection/acceptance, but is not able 
to say anything about the relative probability of certain 
regions of the parameter space with respect to others, 
neither identify a best-fit model given the results ob- 
tained. This is the reason we have chosen to show the 
allowed MSSM space as a uniform shaded area, since in- 
dividual model density in the plot would not carry any 
statistical meaning. Analyses that use experimental re- 
sults to assign probabilities to the SUSY parameter space 
are becoming common. A first evaluation of the possi- 
bilities that IceCube presents in probing the constrained 
MSSM, including individual model information, was done 
in Ref. [41]. A more elaborate analyses in this direc- 
tion, including event energy and direction information, 
as well as constraints from direct dark matter search ex- 
periments and accelerator results in a global fit to the 
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FIG. 6. Left: 90% confidence level upper limits on the spin-dependent neutralino-proton cross section, cr^p , as a function of 
neutralino mass. The results from the analyses presented in this paper are shown as the black dots joined with lines to guide the 
eye (solid and dashed for the W^W~ and the bb annihilation channels respectively). The shaded area represents the allowed 
MSSM parameter space taking into account current accelerator, cosmological and direct dark matter search constrains. The red 
curve shows the expected sensitivity of the completed IceCube detector. Results from Super-K [42] , KIMS [43] and COUPP [44] 
are also shown for comparison. Right: 90% confidence level upper limits on the spin-independent neutralino-proton cross 
section, a^p, as a function of neutralino mass. The results from the analyses presented in this paper are shown as the black 
dots joined with lines to guide the eye (solid and dashed for the W^W~ and the bb annihilation channels respectively). The 
shaded area represents the allowed MSSM parameter space taking into account current accelerator, cosmological and direct 
dark matter search constrains. The red curve shows the expected sensitivity of the completed IceCube detector. Results from 
CDMS [36] and Xenon [37] are also shown for comparison. 



SUSY parameter space is being developed by the collab- 
oration with data obtained with the 79-string detector. 

Given that the Sun is essentially a proton target and 
that the muon flux at the detector can be related to the 
capture rate of neutralinos in the Sun, the IceCube lim- 
its on the spin-dependent neutralino-proton cross section 
are currently well below the reach of direct search exper- 
iments, proving that neutrino telescopes are competitive 
in this aspect. For the spin-independent limits, however, 
direct dark matter search experiments can be competitive 
due to the choice of target. Indeed, the latest results from 
the XENONIOO collaboration [37], using 100 days of live- 
time, have already produced stronger spin independent 
limits than those we present in this paper, as shown in 
the right plot of figure 6. However there is some comple- 
mentarity between direct and indirect searches for dark 
matter given the astrophysical assumptions inherent to 
the calculations. Both methods are sensitive to opposite 
extremes of the velocity distribution of dark matter par- 
ticles in the Galaxy (low-velocity particles are captured 
more efficiently in the Sun, high-velocity particles leave 
clearer signals in direct detection experiments), as well 
as presenting different sensitivity to the structure of the 
dark matter halo (a local void or clump can deplete or 
enhance the possibilities for direct detection). 

The data set used in Analysis B covered the time un- 
til the decommissioning of the AMANDA-II detector 
in 2009. The denser configuration of the AMANDA-II 
strings was of key importance on increasing the sensitiv- 
ity to low neutralino masses, while the sparsely spaced 



IceCube strings alone would have yielded a worse result. 
In order to supplant the role of AMANDA-II as a low- 
energy array, the IceCube collaboration has deployed the 
DeepCore array [45] in the clear South Pole ice, in the 
middle of the IceCube layout. DeepCore lies about 500 m 
deeper than AMANDA-II and its placement in the cen- 
ter of IceCube means that three layers of IceCube strings 
can be used as a veto to reject down-going atmospheric 
muons. The deployment of DeepCore was finalized in 
December 2010 and it is currently taking data embedded 
in the IceCube data acquisition system. DeepCore is ex- 
pected to lower the energy threshold of IceCube to the 
0(10 GeV) region, and therefore be an important asset 
in future dark matter searches with IceCube. 
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